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(57) Abstract: Theinvention relates to a method of recon- 
structing a high-resolution 3D image (G) of an examina- 
tion zone of a patient from a 3D image data set (D) of the 
examination zone, the examination zone being subject to 
a periodic motion which is measured, in parallel with the 
acquisition of the 3D image data set (D), as a motion signal 
(E) which represents the periodic motion. In order to en- 
able the formation of high -re solution 3D images having an 
enhanced image quality, the invention proposes a method 
of this kind which includes the following steps: a) recon- 
structing a number of low-resolution 3D images (I) from 
the 3D image data set (D), the low-resolution 3D images 
(I) being reconstructed from 3D image data of the 3D im- 
age data set (D) which has been acquired during different 
phases of motion of the periodic motion, b) determining 
motion information (B) of at least one sub-zone (A) of the 
examination zone during the different phases of motion by 
means of the low-resolution 3D images (I), c) selecting a 
temporal reconstruction window (T) in which the motion 
of the at least one sub-zone (A) is below a predetermined 
level, d) reconstructing a high-resolution sub-image (K) 
for the at least one sub-zone (A) from 3D image data (D) 
lying in the temporal reconstruction window (T) selected 
for the sub-zone (A), and e) forming the desired 3D image 
(G) from the at least one sub-image (K), the zones of the 
3D image (G) which are not reconstructed as sub -images 
then being reconstructed from the 3D image data so as to 
be combined with the at least one sub-image (K). 
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Method of reconstructing a high-resolution 3D image 



The invention relates to a method and a device for the reconstruction of a 
high-resolution 3D image of an examination zone of a patient from a 3D image data set of the 
examination zone, the examination zone being subject to a periodic motion which is 
measured, in parallel with the acquisition of the 3D image data set, as a motion signal which 
represents the periodic motion. 

A method and a device of this kind are known. The reconstruction of 3D 
images of moving objects, for example, of the heart of a patient, is nowadays carried out by 
selecting a suitable temporal reconstruction window in which acceptable results can be 
achieved for the overall anatomy to be imaged. This means that from the 3D image data set 
acquired that 3D image data which has been acquired in a given time window in which, for 
example, the anatomy to be imaged has been subject to the least motion, is selected for the 
reconstruction. In the case of aperiodic motion of the anatomy, for example, the cardiac 
motion or the respiratory motion, therefore, only 3D image data from a given phase of 
motion are used for the reconstruction, whereas all other reconstruction data acquired is not 
evaluated. For the imaging of the heart or the coronary vessels, for example, it is preferred to 
use exclusively 3D image data acquired during the diastole, whereas the formation of 3D 
images of the abdomen utilizes exclusively 3D data which has been acquired in the state of 
exhalation. 

The motion of given parts of organs or notably of different parts of the heart, 
however, does not take place simultaneously but follows the excitation pattern of the organ; 
this means that individual parts of the organ to be imaged can move to a different extent at 
different instants. The use of a fixed (temporal) reconstruction window for the entire object to 
be imaged, therefore, often gives rise either to a degraded image quality or to sub-optimum 
use of the data acquired during the assumed phases of rest of the organ. 

Therefore, it is an object of the invention to enhance the method and the 
device of the kind set forth in such a manner that high-resolution 3D images of a moving 
examination zone of a patient can be formed with an enhanced image quality. 
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This object is achieved in accordance with the invention by means of a method 
as disclosed in claim 1 and by means of a device as disclosed in claim 6. The method is 
characterized in that it includes the following steps: 

a) reconstructing a number of low-resolution 3D images from the 3D image data 
5 set, the low-resolution 3D images being reconstructed from 3D image data of the 3D image 

data set which has been acquired during different phases of motion of the periodic motion, 

b) determining motion information of at least one sub-zone of the examination 
zone during the different phases of motion by means of the low-resolution 3D images, 

c) selecting a temporal reconstruction window in which the motion of the at least 
10 one sub-zone is below a predetermined level, 

d) reconstructing a high-resolution sub-image for the at least one sub-zone from 
the 3D image data lying in the temporal reconstruction window selected for the sub-zone, and 

e) forming the desired 3D image from the at least one sub-image, the zones of the 
3D image which are not reconstructed as sub-images then being reconstructed from the 3D 

15 image data so as to be combined with the at least one sub-image. 

Claim 6 discloses the corresponding device which is provided with a 
reconstruction unit as well as an arithmetic unit in order to carry out the method. 

The invention is based on the recognition of the fact that for the reconstruction 
of a 3D image it is less than optimum to use only 3D image data from the same temporal 

20 zone, that is, from a fixed phase of motion, for the entire examination zone, for example, an 
organ to be examined such as the heart Because individual parts of the examination zone 
could move to a different extent at different instants or in different phases of motion, in 
accordance with the invention it is proposed to select specific reconstruction windows for the 
individual parts of the examination zone which move differently, to reconstruct sub-images 

25 in such reconstruction windows, and to combine said sub-images subsequently so as to form 
an overall image, that is, the desired 3D image. Each individual sub-image can thus be 
derived from different 3D image data of the same 3D image set. 

The complete examination zone to be reproduced in a 3D image can thus be 
subdivided into individual sub-zones for which a respective optimum (temporal) 

30 reconstruction window is determined each time in accordance with the invention. However, 
individual sub-zones of the examination zone which are of particular interest can also be 
selected in order to perform such a formation of sub-images therein on the basis of an 
adapted reconstruction window, whereas all other zones of the examination zone are 
reconstructed in the conventional manner on the basis of all 3D image data of the 3D image 
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data set or on the basis of a special sub-quantity thereof which was acquired in a low-motion 
phase. 

Further versions and preferred embodiments of the method and the device in 
accordance with the invention are disclosed in the claims 2 to 5 and 7 and 8, respectively. 
5 The invention is particularly suitable for use in the reconstruction of 3D images of the heart 
as well as of the coronary vessels; in that case the motion signal representing the periodic 
cardiac motion preferably corresponds to an electrocardiogram. The invention, however, in 
principle can also be used for the imaging of other objects or other zones, for example, for 
the imaging of the abdomen or of individual organs such as the liver or the stomach, because 

10 such zones or organs are also subject to a periodic motion, that is, the respiratory motion. A 
respiratory motion signal can then be measured as the motion signal; various methods are 
suitable for this purpose, for example, a signal which represents the motion of the abdominal 
wall or of the diaphragm. Carrying out the invention in principle offers an enhanced image 
quality, notably when individual sub-zones or sub-objects of the examination zone to be 

1 5 imaged exhibit a non-simultaneous pattern of motion. 

According to a preferred version of the invention, a motion model of the 
examination zone is formed in order to determine the motion of the at least one sub-zone; to 
this end, notably significant points in the low-resolution 3D images are tracked or objects 
reproduced in the low-resolution 3D images are segmented. Significant points may be, for 

20 example, anatomical features such as, for example, bifurcations in the coronary vessels. 
Tracking such significant points in time in the individual low-resolution sub-images thus 
enables the acquisition of motion information on how said points have moved in time. 
Alternatively, this motion information can be automatically extracted by way of 
segmentation, for example, of the vascular tree in the low-resolution 3D images associated 

25 with various phases of motion. 

According to a further preferred version a separate temporal reconstruction 
window is selected for each volume element of the examination zone or for regions of the 
examination zone which are of particular interest, each volume element (voxel) or the regions 
of special interest being separately reconstructed on the basis of the 3D image data acquired 

30 in the selected temporal reconstruction window. In as far as adequate calculation time is 

available, or if a 3D image with the highest possible resolution and the best possible motion 
correction is to be formed, the method in accordance with the invention can thus be applied 
to every individual voxel of the examination zone; this means that for each individual voxel 
there is derived motion information on the basis of which a special temporal reconstruction 
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window is determined for the voxel, that is, a window in which the reconstruction of this 
voxel image value takes place. However, it is also possible to combine each time several 
voxels so as to form individual groups of voxels for which the method in accordance with the 
invention is carried out so as to save calculation time. 
5 The invention can be used in principle for the formation of a 3D image from 

arbitrary 3D image data; this means that in principle it is not important which modality is 
used to acquire the 3D image data. The method in accordance with the invention, however, is 
preferably used in conjunction with 3D image data acquired by means of computed 
tomography or by means of a 3D rotation X-ray method, for example, by means of a C-arm 
10 X-ray system. 

The invention also relates to a computer program as disclosed in claim 9 for 
the execution of the method of claim 1 and/or for controlling a device as claimed in claim 6. 
This computer program includes notably programming means which are suitable for carrying 
out the method in accordance with the invention or for controlling the device in accordance 
1 5 with the invention when the computer program is executed by a computer or a suitable 
arithmetic unit. 

The invention will be described in detail hereinafter with reference to the 
drawings. Therein: 

Fig. 1 is a diagrammatic representation illustrating the method in accordance . 
20 with the invention, and 

Fig. 2 shows a block diagram illustrating a device in accordance with the 

invention. 

The method in accordance with the invention will be described in detail 
hereinafter on the basis of an example where a high-resolution 3D image of the vascular tree 

25 of the coronary vessels is to be reconstructed. When the contraction pattern of the heart is 
studied, it appears that this pattern commences with the contraction of the atria, followed by 
the contraction of the ventricles, starting with the apex. Because of the complex motion 
pattern of the anatomy observed, the various volume elements (voxels) along the coronary 
artery exhibit significantly different patterns of motion during the cardiac rhythm or a cardiac 

30 motion phase. Whereas in the known reconstruction method the temporal reconstruction 
window is defined in a simple manner by a constant temporal delay in respect of the R 
deflection in the electrocardiogram (ECG) and a given temporal length of the reconstruction 
window. The same reconstruction window is then used for the reconstruction of all voxels in 
the volume to be reconstructed. 



WO 02/103639 



PCT/IB02/02267 



5 

In accordance with the invention, however, voxel-specific reconstruction 
windows are to be selected; these windows are used first for the reconstruction of individual 
sub-images which are subsequently combined so as to form a desired overall image. Fig. 1 
shows an electrocardiogram E which shows the characteristic variation as a function of time t 
5 of a motion signal representing the cardiac motion (a voltage is measured). A significant 
aspect of the ECG is the so-called R deflection which will be used hereinafter so as to define 
the reconstruction windows. 

Using an imaging system, for example, a computed tomography apparatus or a 
C-arm X-ray device operating on the principle of 3D rotation angiography, first a 3D image 

10 data set D is acquired. This means that individual data set elements, that is, so-called 

projections Doi, D02, of the examination zone are acquired from different perspectives or 
imaging positions, that is, continuously in time t and at equidistant instants. These individual 
image data set elements together form a 3D image data set D wherefrom a 3D image of the 
examination zone can be formed. Because the acquisition of the 3D image data D takes place 

15 temporally in parallel with the acquisition of the ECG E, the 3D image data D acquired can 
be assigned to individual phases of motion of the examination zone. The data set elements 
Doi, Dn, D21, D31, ... have thus all been acquired in a first phase of motion in which the ECG 
E exhibits the R deflection; the data set elements D02, D i2 , D22, D32, ... have been acquired in 
a second phase of motion with a first, fixed temporal delay with respect to the R deflection, 

20 etc. The R deflections occur at a distance in time (period) amounting to T. 

In accordance with the invention respective low-resolution 3D images I are 
formed from the data set elements acquired each time during the same phase of motion. This 
means that, for example, a first low-resolution image Ii is reconstructed from the data set 
elements D 0 i, Dn, D21, D 3 i, ... of the first phase of motion, that a second low-resolution 

25 image I 2 is reconstructed from the data set elements D02, Di 2 , T>22> D 32 , ... of the second phase 
of motion, etc. Each of these low-resolution images l\ 9 I 2 , 13, ... is thus formed from 3D image 
data exhibiting a fixed temporal delay relative to the R deflection. 

From such low-resolution 3D images Ii, I 2 , 13, ... there is derived motion 
information B which provides information on the motion of the anatomy present in the 

30 examination zone during the cycle of motion. In the example described herein, therefore, 
motion information B is to be acquired as to how individual zones of the heart or individual 
coronary vessels move during the cardiac cycle. To this end, for example, significant points 
such as bifurcations of the coronary vessels can be tracked in the individual low-resolution 
images I or the vascular tree can be automatically segmented in these images. It is thus 
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possible to derive infoimation defining the phases of motion in which sub-zones of the 
examination zone, or in an extreme case individual voxels of the examination zone, have 
moved and to what extent and also the phases of motion or periods of time in which no or 
only insignificant motion has occurred during a motion cycle. 
5 This motion information is subsequently used to determine a respective 

optimum temporal reconstruction window Ti, T2, T3, ... for individual sub-zones Ai, A2, A3, 
which may also correspond to individual voxels in an extreme case. For a first sub-zone 
Ai a first temporal reconstruction window Ti is thus determined; this means at the same time 
that only image data set elements acquired from this segment of time of the motion cycle are 

10 used for the reconstruction of a high-resolution sub-image Ki of this sub-zone Ai. This 

reconstruction window Ti can also be recognized in Fig. 2; it can be deduced therefrom that, 
for example, the reconstruction of a high-resolution sub-image Ki of the sub-zone Ai utilizes 
only the image data set elements D02, D03, D04, D 0 s, Dj 2 , D13, Dh, D15, D22, D23, ... because r 
this sub-zone obviously has not been subject to motion, or only hardly so, during the time 

15 window Ti of the motion cycle. For a second sub-zone A2, however, a different temporal 
reconstruction window T2 is derived from the motion information B. For the formation of a 
high-resolution image K 2 of the sub-zone A2, therefore, image data set elements D04, Dos, 
D06, Do7 3 D14, Di5, ... are used, because obviously the sub-zone A2 has moved the least during 
the segment of time T 2 . 

20 For each sub-zone A for which the method in accordance with the invention is 

to be used, therefore, an individual temporal reconstruction window T is selected, the 
temporal delay relative to the R deflection and the duration in time of the reconstruction 
window may be different for each sub-zone A. The reconstruction windows, however, may 
overlap or also be identical. 

25 As has already been described, in an extreme case a separate temporal 

reconstruction window T can be determined for each individual voxel of the examination 
zone, so that each individual voxel is optimally reconstructed individually from specially 
selected data set elements. In another extreme case a separate reconstruction window may be 
determined for a single sub-zone only of the examination zone, whereas all other zones of the 

30 examination zone are reconstructed in a customary manner from all 3D image data of the 3D 
image data set D or from 3D image data of a fixed temporal reconstruction window. It is to 
be noted again that the reconstruction windows are variable and cover for each sub-zone a 
different length and a different segment in time of the motion period. 
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In a final step the high-resolution sub-images Ki, K2, K3, ... formed are 
combined so as to form the desired high-resolution 3D image G. The sub-zones in this 3D 
image for which a separate reconstruction window was selected in accordance with the 
invention thus exhibit a clearly enhanced image quality with a higher resolution and fewer 
5 artifacts- 
Fig. 2 is a representation in the form of a block diagram of a device in 
accordance with the invention for carrying out the described method. Therein, first a 3D 
image data set D is acquired by means of an image data acquisition device 1, for example, a 
computed tomography apparatus or a C-arm X-ray device. In parallel therewith a motion 
1 0 signal E, for example, an ECG is acquired by means of a suitable measuring device 2, for 
example, by means of an electrocardiograph; this motion signal provides information as 
regards the motion of the anatomy to be imaged. The data acquired is first stored in a 
memory 3 and subsequently applied to the reconstruction unit 4 for evaluation and further 
processing. 

1 5 Therein, the low-resolution 3D images I, representing the examination zone in 

different phases of motion, are formed first in a first reconstruction module 41 . Subsequently, 
said low-resolution images are applied to an arithmetic unit 5 in which the motion 
information B, for example, in the form of a motion model, is acquired in a first module 5 1 ; 
this motion information is subsequently evaluated in a second module 52. In the module 52, 

20 therefore, the motion information is used to determine the individual temporal reconstruction 
windows T for individual sub-zones A which are defined, for example by a user. These 
reconstruction windows are then applied to a second reconstruction module 42 which forms 
the respective high-resolution sub-images K for the individual sub-zones A on the basis of 
this information and as well as on the basis of the motion signal E from the 3D image data 

25 set. These high-resolution sub-images K are ultimately combined in a combination module 
43 so as to form the desired high-resolution 3D overall image G of the examination zone. 
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CLAIMS: 



1 . A method of reconstructing a high-resolution 3D image (G) of an examination 
zone of a patient from a 3D image data set (D) of the examination zone, the examination zone 
being subject to a periodic motion which is measured, in parallel with the acquisition of the 
3D image data set (D), as a motion signal (E) which represents the periodic motion, which 

5 method includes the steps of: 

a) reconstructing a number of low-resolution 3D images (I) from the 3D image 

data set (D), the low-resolution 3D images (I) being reconstructed from 3D image data of the 
3D image data set (D) which has been acquired during different phases of motion of the 
periodic motion, 

10 b) determining motion information (B) of at least one sub-zone (A) of the 

examination zone during the different phases of motion by means of the low-resolution 3D 
images (I), 

c) selecting a temporal reconstruction window (T) in which the motion of the at 

least one sub-zone (A) is below a predetermined level, 
15 d) reconstructing a high-resolution sub-image (K) for the at least one sub-zone 

(A) from 3D image data (D) lying in the temporal reconstruction window (T) selected for the 
sub-zone (A), and 

e) forming the desired 3D image (G) from the at least one sub-image (K), the 

zones of the 3D image (G) which are not reconstructed as sub-images then being 
20 reconstructed from the 3D image data so as to be combined with the at least one sub-image 
(K). 

2. A method as claimed in claim 1, characterized in that a motion model of the 
examination zone is formed in order to determine the motion information (B) of the at least 

25 one sub-zone (A), that is, notably by tracking significant points in the low-resolution 3D 
images or by segmentation of objects reproduced. 



3. A method as claimed in claim 1, characterized in that for each voxel of the 

examination zone, or for regions of the examination zone which are of particular interest, 
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there is selected a separate temporal reconstruction window (T), and that each volume 
element or each region of particular interest, is separately reconstructed on the basis of the 
3D image data acquired in the selected temporal reconstruction window (T). 

5 4. A method as claimed in claim 1, characterized in that the 3D image data set 

(D) has been acquired by means of a computed tomography apparatus or a 3D rotation X-ray 
device. 

5. A method as claimed in claim 1, characterized in that an electrocardiogram or 
10 a respiratory motion signal is used as the motion signal (E) representing the motion of the 

examination zone. 

6. A device for reconstructing a high-resolution 3D image (G) of an examination 
zone of a patient from a 3D image data set (D) of the examination zone, the examination zone 

15 being subject to a periodic motion which is measured, in parallel with the acquisition of the 
3D image data set (D), as a motion signal (E) which represents the periodic motion, which 
device includes 

a) a reconstruction unit (4) for the reconstruction of a number of low-resolution 
3D images (I) from the 3D image data set (D), the low-resolution 3D images (I) being 

20 reconstructed from 3D image data of the 3D image data set (D) which has been acquired 
during different phases of motion of the periodic motion, 

b) an arithmetic unit (5) for determining the motion of at least one sub-zone (A) 
of the examination zone during the different phases of motion by means of the low-resolution 
3D images (T), 

25 the arithmetic unit being arranged to select a temporal reconstruction window (T) in which 
the motion of the at least one sub-zone (A) is below a predetermined level, and 
the reconstruction unit (4) being arranged to reconstruct at least one high-resolution sub- 
image (K) for the at least one sub-zone (A) from 3D image data (D) lying in the temporal 
reconstruction window (T) selected for the sub-zone (A), and to form the desired 3D image 

30 (G) from the at least one sub-image (K), the zones of the 3D image (G) which are not 

reconstructed as sub-images then being reconstructed from the 3D image data so as to be 
combined with the at least one sub-image (K). 
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7. A device for forming 3D images of an examination zone of a patient which 
includes an image acquisition device (1) for the acquisition of a 3D image data set of the 
examination zone, the examination zone being subject to a periodic motion, and also includes 
a measuring device (2) for measuring, in parallel with the acquisition of the 3D image data 

5 set, a motion signal which represents the periodic motion, and also a device for the 
reconstruction of a high-resolution 3D image as claimed in claim 6. 

8. A device as claimed in claim 7, characterized in that the device is a computed 
tomography apparatus or a 3D rotation X-ray device, notably a C-arm X-ray device. 

10 

9. A computer program for carrying out the method claimed in claim 1 and/or for 
controlling a device as claimed in claim 6. 
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